Constraining dark matter annihilation with the isotropic 7-ray background: 

updated limits and future potential 



Torstcn Bringmanro and Francesca Calorqt| 
II. Institute for Theoretical Physics, University of Hamburg, 
Luruper Chaussee 149, 22761 Hamburg, Germany 

Mattia Di Maurc@ and Fiorenza Donate! 
Dipartimento di Fisica, Universita di Torino and IN FN, 
Sezione di Torino, Via P. Giuria 1, 10125 Torino, Italy 

The nature of the Isotropic 7-ray Background (IGRB) measured by the Large Area Telescope 
(LAT) on the Fermi 7-ray space Telescope {Fermi) remains partially unexplained. Non-negligible 
contributions may originate from extragalactic populations of unresolved sources such as blazars, 
star-forming galaxies or galactic milli-second pulsars. A recent prediction of the diffuse 7-ray emis- 
sion from Active Galactic Nuclei (AGN) with a large viewing angle with respect to the line-of-sight 
(l.o.s.) has demonstrated that this faint but numerous population is also expected to contribute 
significantly to the total IGRB intensity. A more exotic contribution to the IGRB invokes the pair 
annihilation of dark matter (DM) weakly interacting massive particles (WIMPs) into 7-rays. In 
this work, we evaluate the room left for galactic DM at high latitudes (> 10°) by including photons 
from both prompt emission and inverse Compton scattering, emphasizing the impact of the newly 
discovered contribution from misaligned AGN (MAGN) for such an analysis. Summing up all signif- 
icant galactic and extragalactic components of the IGRB, we find that an improved understanding 
of the associated astrophysical uncertainties is still mandatory to put stringent bounds on thermally 
produced DM. On the other hand, we also demonstrate that the IGRB has the potential to be one of 
the most competitive future ways to test the DM WIMP hypothesis, once the present uncertainties 
are even slightly reduced. In fact, if MAGN contribute even at 90% of the maximal level consistent 
with our current understanding, thermally produced WIMPs would be severely constrained as DM 
candidates for masses up to several TeV. 

PACS numbers: 95.30.Cq,95.35+d,95.85.Pw,96.50.sb 



I. INTRODUCTION 

The nature of dark matter (DM) remains one of the 
most intriguing mysteries in fundamental physics and 
cosmology. A viable investigation technique consists in 
the search for its stable annihilation products in the halo 
of galaxies, and in particular in the Milky Way. Mak- 
ing the common hypothesis for DM to consist of weakly 
interacting massive particles (WIMPs), one of the most 
promising indirect detection channels is through annihi- 
lation into 7-rays (see Ref. [l| for a recent review). Un- 
precedented 7-ray observations of diffuse emission and 
point sources by the Fermz-LAT [2j have stimulated the 
search for exotic components from DM annihilation in 
the Milky Way. in extragalactic nearby objects, as well 
as in cosmological structures [H-Q- The incontrovert- 
ible DM signature would be a 7-ray monochromatic line 
at energies around the WIMP mass. A recent analysis 
of Fermz'-LAT data of the galactic center region has un- 
veiled the feature of a narrow line close to 130 GeV [HHI, 
with an angular distribution consistent with the expecta- 
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tion for WIMP DM [l[ and a statistical significance of at 
least 3tr after trials factor correction (in fact, more than 
5(7 when taking into account both spectral and spatial 
information of the incoming events |ll| ) - though results 
from the Fermz'-LAT Collaboration investigating this ef- 
fect, albeit in somewhat different target regions, are less 
optimistic [l2l Il3| . The interpretation in terms of annihi- 
lating DM would in any case require the loop-suppressed 
direct annihilation into photons to be larger than typi- 
cally expected for thermally produced DM, a fact which 
has already inspired considerable model-building efforts 
(see, e.g., the compilation in Refs. [HtHI 101 an overview). 
Only a larger exposure will allow to finally settle possible 
issues concerning the statistical robustness and potential 
instrumental effects of the line feature; with the proposed 
new observation strategy for Fermi, this could be achiev- 
able until the end of 2014 

A faint 7-ray emission at high latitudes (|6| > 10°) has 
been measured by Fermz'-LAT 15|, showing a high de- 
gree of isotropy [lfj and an energy spectrum compatible 
with a power-law. Though its very nature is still unex- 
plained, it is reasonable to believe it to be the superpo- 
sition of (mainly) extragalactic unresolved point-sources 
and, to a lesser extent, of astrophysical diffuse processes 
(see Ref. [ljj and Section [TT| for a detailed discussion). A 
possible contribution from DM annihilation, both from 
galactic and cosmological distances, can be hidden in this 
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isotropic 7-ray background (IGRB) data as well fl7l - l25l | . 1 
Previous conservative upper limits on the DM annihila- 
tion cross section from the high latitudes 7-ray emission 
were derived e.g. in [T^], by subtracting from the IGRB 
[l5| the minimum estimated fluxes for the most signifi- 
cant unresolved source populations. 
The faint but numerous population of Active Galactic 
Nuclei (AGN) whose jet axes have a large viewing angle 
with respect to the line-of-sight (l.o.s.) has recently been 
shown [27], [H] to contribute between 10% and 100% to 
the IGRB intensity, depending on the astrophysical un- 
certainties involved. In the present paper, we extend and 
update the analysis in Ref. [l7j by taking into account 
those new results on the isotropic diffuse flux predicted 
for misaligned AGN (MAGN) ^EH- We find that the 
annihilation intensity of DM in the galactic halo may be 
severely constrained by the inclusion of this new diffuse 
astrophysical background, depending on the WIMP mass 
and on the annihilation channel considered. 

In this context, it is worth noting that classical DM in- 
direct detection targets like the galactic center or dwarf 
spheroidal galaxies, which currently result in the best 
available limits on DM annihilation 0, Hi| , will at some 
point inevitably face fundamental limitations in the sense 
that even greatly increased statistics will not allow to 
further improve those limits [l|. For ground- based tele- 
scopes, this is likely the case already for anything beyond 
the planned Cherenkov Telescope Array [3fJ , while for fu- 
ture space-based missions there might be room for further 
improvement of the limits by up to one order of magni- 
tude with respect to what can be expected by the end of 
the Fermi mission. Once the MAGN-related astrophysi- 
cal uncertainties are under better control, the IGRB may 
thus well turn out to be the new driving force in setting 
ever more stringent limits on the DM annihilation rate. 

The paper is organized as follows. In Scct.[TT]we present 
the possible contributions to the IGRB and we discuss 
the relevance of unresolved extragalactic sources, notably 
MAGN. The expected galactic DM signal is described in 
Sect. IIII1 The DM spectrum is made up by the prompt 7- 
ray component and photons from inverse Compton scat- 
tering (ICS) of electrons and positrons produced by DM 
annihilation. By subtracting the astrophysical compo- 
nents from the high latitudes 7-ray emission, in Scct. lIVI 
we derive upper limits on the DM annihilation cross sec- 
tion for several annihilation channels and under various 
assumptions for the MAGN contribution. We discuss 
the implications for DM indirect detection, as well as the 
dependence on the other possible astrophysical contribu- 
tions. In Sect. [V] we present our conclusions. 



The 130 GeV feature described above does not show up in 
the IGRB — as expected for the annihilation rate of (cru)-^ ~ 
10~ 27 cm 3 s _1 inferred from the galactic center observation Q 
(unless one adopts very optimistic models for the evolution and 
distribution of DM subhalos [26l0. 



II. DIFFUSE 7-RAY EMISSION FROM 
UNRESOLVED ASTROPHYSICAL SOURCES 

The IGRB found at high latitudes (|6| > 10°) has 
been obtained by the Fermi-LAT Collaboration |15| af- 
ter subtracting from the LAT data the contributions 
from point-like sources, the galactic diffuse emission 
(which represents the major systematic uncertainty of 
that analysis), the background from CRs in the detec- 
tor and solar photons. The energy spectrum of the 
IGRB is consistent with a power-law having spectral in- 
dex a = 2.41 ± 0.05 and integrated intensity /(> 100 
McV) = (1.03 ±0.17) • 10~ 5 cm" 2 s" 1 sr" 1 . This diffuse 
emission likely arises from the superposition of several 
contributions from unresolved point-source populations, 
as well as from photons originating from truly diffuse pro- 
cesses (i.e extended source contributions like from ultra- 
high energy cosmic rays or gravitationally induced shock 
waves during structure formation, see Ref. (l7| and ref- 
erences therein). Most of the photons in the IGRB are 
believed to be of extragalactic origin. The truly diffuse 
processes are subdominant and may be non-negligible (% 
level) only at the high-energy end of the range explored 
by the Fermi-LAT . Given their very marginal role in the 
analysis undertaken here, they will be neglected from now 
on. 

Misaligned A GN. Recently, a thorough analysis of the 7- 
ray flux from unresolved MAGN has demonstrated that 
this numerous but faint population may contribute be- 
tween 10% and 100% to the measured IGRB [H|2|. A 
jet displacement of about 14° marks the separation be- 
tween blazars and non-blazar (namely misaligned) AGN. 
In the unification model [3l[, radio galaxies are those 
misaligned active galaxies whose jet points on average 
at an angle > 44°, while lower angles indicate typically 
steep spectrum radio quasars. The result in Refs. [27U28j 
is based on the 7-ray data of 12 MAGN by Fermi-LAT, 
and is based on a correlation between their 7-ray and 
core radio luminosities. This correlation has been tested 
against the MAGN number count distribution, and val- 
idated by 95% flux upper limits on the Fermi-LAT 7- 
ray data for a large sample of bright core radio MAGN. 
These results support the hypothesis that 7-ray photons 
are produced by synchrotron-self Compton scatterings 
by the same electrons responsible for the synchrotron 
emission at GHz frequencies. Ref. [27j quantifies sev- 
eral sources of uncertainty and estimates them to affect 
the diffuse emission from unresolved MAGN by about a 
factor of 10. The expected flux is therefore provided as 
a band (roughly decreasing with E~ 2 - 35 ) affected by the 
extragalactic background light (EBL) absorption above 
50 GeV. At 1 GeV, the MAGN flux (multiplied by E 2 ) 
ranges from 6-10" 7 GeV cm" 2 s" 1 sr" 1 to 5-10" 6 GeV 
enr 2 s" 1 sr -1 , while at 100 GeV the EBL absorption 
starts to be effective and the flux spans from 9T0 -9 GeV 
cm" 2 s" 1 sr -1 to 9-10 -8 GeV cm" 2 s" 1 sr" 1 . 
The result in Ref. {27} states that MAGN account for 
at least 10% of the IGRB measured by Fermi-LAT [H| 
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and can in principle explain most of it. Added to other 
possible diffuse emission components, this might almost 
close the room to an exotic diffuse contribution such as 
from DM annihilation in space. This point is the focus of 
our paper: evaluating how much the MAGN 7-ray diffuse 
emission constrains the presence of WIMP DM annihilat- 
ing in the galactic halo. In the remainder of this Section, 
we briefly discuss the astrophysical sources leading to a 
non-negligible diffuse 7-ray emission at \b\ > 10° (see 
Ref. 17J for further details). 

Blazars. The census of the sources detected by the Fermi- 
LAT [32l I331] indicates that the largest population is made 
by blazars, a class of AGN with their jet oriented along 
the l.o. s., thus leading to a significant observed 7-ray flux. 
The analysis of the spectral and statistical properties of 
the resolved Fermi-LAT blazars [34| shows that the un- 
resolved counterpart is a sizable component, able to ex- 
plain up to 20% - 30% of the IGRB. The contribution 
from the two distinct populations of Flat Spectrum Ra- 
dio Quasars (FSRQs) and BL Lacertae (BL Lacs) objects 
has been found to follow a power-law spectrum with re- 
spectively a somewhat softer (o-fshQs = 2.45 ±0.03) and 
harder (aBLLacs = 2.23 ± 0.03) spectral index. Recently, 
Ref. [35| has evaluated the FSRQ contribution to the dif- 
fuse 7-ray flux using a new determination of the spectral 
energy distribution and luminosity function. The result- 
ing flux is predicted with an uncertainty of less than a 
factor of two. It represents no more than 10% of the 
IGRB in the 0.1-100 GcV, and becomes negligible for E 
> 20 GcV. 

Milli-Second Pulsars. The second-most abundant pop- 
ulation of the Fermi-LAT sky is represented by galactic 
pulsars and Milli-Second Pulsars (MSPs). It has been 
argued that the latter population may be brighter in 7- 
rays than ordinary pulsars as a consequence of their rapid 
spin frequency [3g|. Moreover, MSPs are expected to be 
more numerous at high latitudes, since they have typical 
ages that can largely exceed the oscillation time across 
the galactic disk [37| ■ The prediction of the diffuse 7-ray 
flux due to high latitudes unresolved MSPs is estimated 
to be about 10% at ~ 3 GeV, where it reaches its max- 
imum value due to the adopted broken power-law spec- 
trum 17]. 2 In the following, we refer to the unresolved 
MSPs diffuse 7-ray flux and its estimated uncertainty as 
reported in Ref. [17] . 

Star-forming galaxies. Though their number density is 
larger than that of the blazars, star-forming (SF) galax- 
ies are more difficult to detect in 7-rays because of their 
less intense luminosities. Based on three years of Fermi- 



2 In the Second Fermi catalog [3j| only about 20 MSPs at |6| > 
10° are listed, but the number of galactic (non-globular cluster) 
MSPs has recently risen to 40 |3SH - The statistical properties of 
such a population could improve the prediction for the unresolved 
component. Such an analysis is beyond the scope of the present 
paper, but is anyway not expected to modify any of the results 
and conclusions presented here. 



LAT data, Ref. [39J analyzes a sample of 64 SF galaxies 
beyond the Local Group, detecting only four objects in 7- 
rays: two are typical star-burst galaxies (M82 and NGC 
253), while the other two host Seyfert 2 nuclei (NGC 
1068, NGC 4945). None of the remaining 60 source can- 
didates has shown an excess above the background. The 
contribution of SF galaxies to the IGRB is evaluated in 
Ref. (39( | by exploiting the correlation between infrared 
and 7-ray luminosities [40j |. In order to compute the dif- 
fuse emission, the differential photon flux of an individual 
galaxy has to be modeled. Given the lack of statistics, 
the 7-ray average spectral properties of SF galaxies are 
difficult to firmly establish, although differences in the 
spectra of quiescent and star-burst galaxies are expected 
to be present. To take into account this effect, Ref. [39| 
adopts two possible 7-ray spectra supposed to bracket the 
possible contribution from multiple types of SF galaxies. 
The first one refers to Milky Way-like SF galaxies (model 
MW), the second one assumes a power- law spectrum in 
agreement with the properties of the Fermi-LAT detected 
star- burst galaxies (model PL). The predictions differ in 
particular above 5 GcV, where the MW model gets sig- 
nificantly softer. At 100 GeV the PL model overprcdicts 
the MW one by one order of magnitude. The estimates 
span 4%-23% of the IGRB intensity above 100 MeV. 
As a summary of this Section, we show in Fig. Q] the con- 
tributions to the IGRB discussed above, with uncertainty 
bands as reported in the relevant references (roughly cor- 
responding to la deviations). The Fermi-LAT IGRB 
data [l5l | are shown along with the power-law fit (solid 
black line). The solid dark-green curve refers to the ex- 
pected diffuse 7-ray flux from unresolved MAGN f27j . 
the light-green band bracketing the corresponding uncer- 
tainty. The other displayed components are: SF galaxies 
(the MW model is depicted by the dashed dark- (light-) 
blue line (area)), MSPs (dot-dashed red line and orange 
band), FSRQs (purple dotted line and pink band), and 
BL Lacs (dotted dark-grey line and light-grey band). 
The BL Lacs flux follows a power-law, while the specific 
shape of the FSRQ contribution is the result of a more re- 
cent, sophisticated analysis of their spectral energy distri- 
bution and luminosity function (3||. The crossing point 
for the two curves is around 1-2 GeV: above this energy 
the BL Lacs flux dominates over the one from FSRQs. 
The 7-rays from unresolved MSPs show a peculiar spec- 
trum peaked at about 1 GcV (resulting from the primary 
electron cooling due to curvature radiation) and domi- 
nate over the blazar spectra from 300 MeV up to 3-4 
GeV. A similar low-energy bumped flux comes from the 
SF galaxies contribution. The MAGN is clearly the domi- 
nant component over the whole energy range (notice that 
for E 1 > 30 GcV the BL Lacs flux is actually expected 
to drop in the same way as the MAGN flux, due to EBL 
absorption, but this has not been taken into account in 
the flux estimate referred to here). In the following, we 
define Sbms as the sum of the contributions from Blazars 
(both BL Lacs and FSRQs), MSPs and SF galaxies. By 
varying the contribution of each component within its 
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FIG. 1: The diffuse 7-ray fluxes predicted for various unre- 
solved point-source populations, along with the IGRB data 
from the Fermi-LAT Collaboration based on high galactic 
latitude > 10°) observations [l^]. The displayed con- 
tributions to the IGRB and corresponding uncertainty bands 
arise from: MAGN [23] (solid dark-green curve and light- 
green band), MW model for SF galaxies [3^| (dot-dashed red 
line and orange band) , MSPs [l7l ] (dot-dashed red line and or- 
ange band) , BL Lac s |34ll (dotted dark-grey line and light-grey 
band), and FSRQs [35| (purple dotted line and pink band). 



uncertainty band, the IGRB data can be saturated by 
multiple combinations of the above fluxes. Alternatively, 
only little room is left to other unknown 7-ray sources 
as, e. g., DM. 



III. DIFFUSE 7-RAY EMISSION FROM DARK 
MATTER IN THE GALACTIC HALO 



The self-annihilation of DM particle pairs in the haloes 
of galaxies may give birth, among other species, to 7- 
rays. The intensity of these energetic photons depends 
on the elementary process at stake and on the spatial dis- 
tribution of DM. The flux $ 7 (i? 7 , ip) of 7-rays produced 
by WIMP pair annihilation in the angular direction ip is 
given by gJIH : 



$ 7 (£; 7 ,V) = -r^r-nr-^i 



47T m\ dE 1 2' 



(1) 



where m x is the WIMP mass 



The factor (av) defines 
the annihilation cross section times the relative velocity, 
averaged over the galactic velocity distribution function, 
while dN 1 /dE 1 is the 7-ray production source spectrum 
per DM annihilation event. The last term I(ip) is the in- 
tegral performed along the l.o.s. of the squared DM den- 
sity distribution: 



p 2 (r(A,V))dA. 



(2) 



Here ip is the angle between the l.o.s. and the direction 
towards the galactic center, defined as a function of the 



galactic latitude b and longitude I (cosip = cos b cos I). 
When comparing with experimental data, Eq. ^ must 
be averaged over the telescope viewing solid angle, Afi: 
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I(ip(b,l))dSl. 



(3) 
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The photon spectrum dN^/dE^ depends on the elemen- 
tary processes ruling the annihilation. For the continuum 
7-ray flux, we consider here i) prompt emission, where 
the photons are found in final-state showers or hadronic 
decays of the annihilation products, and ii) ICS by en- 
ergetic electrons and positrons - produced in the same 
way - off the interstellar radiation field (ISRF) . We have 
calculated the spectrum for the prompt emission of both 
photons and e ± with the Pythia Montecarlo code (ver- 
sion 8.162) Q. 

The ICS has been computed from those e ± interacting 
with ambient photons, following the prescription in [45j |. 
Once produced, energetic electrons may diffuse due to 
stochastic scatterings on galactic magnetic field inho- 
mogeneities, and loose energy through interactions with 
the ISRF. The ISRF is composed by the Cosmic Mi- 
crowave Background (CMB) photons, the infrared radi- 
ation (IR) produced by the absorption and the subse- 
quent re-emission of starlight by the galactic dust, and 
by the starlight (SL) originating from stars of the Galac- 
tic disk. Given the energies at stake, galactic diffusion 
may be safely neglected and only energy losses have been 
included, as usually assumed in the literature [HI . Elec- 
tromagnetic energy losses are treated here in the fully 
relativistic Klein-Nishina regime. The photon density 
distribution has been assumed as an average field with 
different normalizations for different sky regions (the SL 
and IR are concentrated along the galactic disk), follow- 
ing Ref. [H]. Specifically, between 10° and 20° all the 
three fields (CMB, IR, SL) are present, while for lati- 
tudes \b\ > 20° only CMB remains. Energy losses due 
to synchrotron radiation have been taken into account in 
the region of latitude 10° < \b\ < 20°, with a value of 
B = 2 fiG for the galactic magnetic field, and neglected 
above \b\ > 20° given the very low intensity of B at high 
latitudes pTl - IHol ] . Switching off the synchrotron radiation 
would increase the ICS 7-ray flux from DM by less than 
30%, resulting in a decrease of the upper limits on (av) 
(shown in the following Sect. IIV|) by 10%. 

The geometrical factor I(ip) in Eg. ([2]) has been evalu- 
ated adopting the Burkert model [51| for the radial DM 
density distribution: 



p(r) 



1+ & 



1 + ^ 



(4) 



The values of r s and p s have been derived fixing the lo- 
cal DM density p(r = r©) = 0.4 GcV cm" 3 (r =8.33 
kpc) and fitting the halo DM mass M(r < r ) contained 
within a certain radius r . U sing available M(r < r ) 
data for different values of tq [53 - [60l | we find a best fit 
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value of r s = 15 kpc. Numerically, we derive for Jao in 
Eq. ©: 2.4, 1.9, 1.7, 2.0 (in units of GeV 2 cm" 6 kpc) 
for the sky regions 10° < \b\ < 20°, 20° < \b\ < 90°, 
|6| > 60°, and |6| > 10° respectively. While in principle 
limits might somewhat improve when taking into account 
the different values of Iaq for different regions of the sky, 
we will in the following only consider the whole region 
|6| > 10° considered by the Ferrm-LAT analysis [l5[ be- 
cause not all of the astrophysical contributions quoted 
above can be easily rescaled to smaller subregions. We 
also note that when integrating Eq. ^ around the galac- 
tic center, different DM distributions p(r) may lead to 
very different results. However, since we study here only 
high latitude regions, the various models for p(r) give 
values for I(ip) which are very similar, and the uncer- 
tainty due to DM density modeling is rather small. For 
example, we find I An (\b\ > 10°) • GeV~ 2 cm 6 kpc _1 = 2.4 
(2.5, 2.8) for the case of an isothermal sphere (a Navarro- 
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an Einasto |62|) profile, which repre- 
sent other typical assumptions about the DM distribu- 
tion. While the existence of a dark disk [63| may further 
enhance those values slightly, we will not consider this 
possibility in the following. We also do not include any 
contribution from sub-haloes, as predicted in numerical 
simulations of cosmological structure formation [64|, [65| , 
and assume a smooth distribution of DM in the galac- 
tic halo. Were the Milky Way halo to be clumpy, the 
geometrical factor in Eq. ([3]) would be increased by the 
sub-structures. In this case, the upper limits we are going 
to present on (av) would be lowered accordingly. 



IV. CONSTRAINTS ON DARK MATTER 
ANNIHILATION INTO 7-RAYS 

The constraints on DM are set, as usual, in terms of the 
velocity averaged annihilation cross section into 7-rays 
(av) as a function of the WIMP mass m x . All the other 
quantities in Eq. ([1} are kept fixed. The galactic DM 
contribution is calculated following the prescriptions out- 
lined in Sect. IIIII We do not include the diffuse 7-ray flux 
from DM in extragalactic structures [66|, H?} , given the 
large model-dependence on the assumed DM (sub-)halo 
profiles and their evolution. If the 7-rays from DM in 
extragalactic structures or sub-haloes in the Milky Way 
halo were added to our fluxes, the ensuing upper limits 
on (av) would be further decreased. In these regards, the 
results we are going to present are meant to be conser- 
vative, i.e. leading to the highest possible value for the 
excluded (av). 

As discussed in Sect. IIIII the geometrical factor at high 
latitudes is very slightly affected by modeling uncertain- 
ties. The choice we made on the Burkert profile and a 
smooth halo is meant to be conservative and different 
assumptions on the DM profile would lower the upper 
bounds wc find by up to about 30% (in fact, considerably 
more for optimistic assumptions about DM substructure 
[68T - l70j h A further uncertainty factor is the value of the 
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FIG. 2: Conservative upper limits on the velocity averaged 
DM annihilation cross section (av) into e final states, as a 
function of the WIMP mass m x . The blue dashed curve cor- 
responds to the 7-ray spectrum from the prompt component 
only, the solid red line refers to the inclusion of ICS. The grey 
area indicates the typically adopted value for the cross section 
of thermally produced WIMPs, (av) ~ 3 ■ 10" 26 cm 3 s" 1 . 



local DM density, which can span the (generous) range 
0.2 - 0.7 GeV cm~ 3 (7lJ. Since p Q enters squared in the 
flux, it has a non-negligible impact on the (av) bounds. If 
not explicitly stated differently, we will take p Q =0.4 GeV 
cm -3 in the following. To keep the analysis as model 
independent as possible, wc will show our results assum- 
ing that DM annihilation proceeds via a single channel 
(branching ratio equal to one), studying separately the 
effect of each final state leading to 7-ray production. 

Finally, the limits on (av) are obtained by requiring 
that the sum of the astrophysical and DM fluxes do not 
exceed any of the Fer?7M-LAT data points fl5j by more 
than 2a, which is an often adopted practice for this kind 
of constraints. We mostly follow this convention, not the 
least in order to allow for a better comparison with previ- 
ous limits, but note that a statistical downward fluctua- 
tions even in a single data bin may lead to a considerable 
underestimation of the allowed DM-induced 7-ray flux if 
the expected astrophysical contributions are very close 
to the observed flux - which, unlike the situation consid- 
ered before in the literature, may indeed be the case when 
taking into account the effect of MAGN. For the case of 
a large MAGN contribution, in particular, wc will there- 
fore also consider the more conservative limits resulting 
from the requirement that no data point is exceeded by 
more than 3a. As can be anticipated from Fig. [TJ those 
limits are always set by the three data bins at 6.4-12.8, 
12.8-25.6 and 25.6-51.2 GeV. 

Fig. [5] represents the exclusion plot for DM annihila- 
tion into e ± and it is emblematic of the effect of the ICS 
inclusion in the DM spectrum. The blue dashed curve 
corresponds to upper limits on (av) when the DM in- 
duced spectrum originates only from the prompt 7-ray 
component, while the solid red line refers to the inclusion 
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FIG. 3: Upper limits on the velocity averaged DM annihilation cross section (an) as a function of the WIMP mass, for different 
levels of Ebms and MAGN. The parameters in the colored areas are allowed by 7-ray data and astrophysical backgrounds by 
the Fermi-ljAT IGRB data [l5|] . In the left panel, we fix Ebms at the minimal level compatible with the expectations and vary 
the contribution from MAGN from the minimal to 70% of the maximally expected level. The right panel shows the impact of 
various choices for Ebms on the limits, fixing the MAGN contribution to the best-fit value derived in Ref. [27j|. The constant 
line at 3 • 10~ 26 cm 3 s -1 represents the 'thermal' cross section of DM in the simplest WIMP models. 



of ICS. The astrophysical background is set by the sum 
of the minimal values for all the components described 
in Sect.|lIJ The effect of ICS inclusion starts to be impor- 
tant for high energies, i. e. at high DM masses, and the 
deviation from the prompt only limits reaches about two 
orders of magnitude at ~ 10 TeV. We notice that the 
uncertainties from ICS are small at the energies at which 
it is effective. For comparison, we also indicate in the 
same plot the 'thermal' cross section of (<jv) ~ 3 ■ 10 -26 
cm 3 s" 1 that results in the correct DM relic density in 
the simplest WIMP scenarios (i.e. velocity-independent 
annihilation and no co-annihilations with other particles 
close in mass to the WIMP). 

In Fig. |21 we study the effect of the various astrophysi- 
cal backgrounds described in Sect.|IT]on the room left for 
a possible DM contribution to the diffuse 7-ray flux. The 
annihilation channel is fixed to bb. In the left panel, we 
study the impact of the previously neglected MAGN com- 
ponent by fixing Ebms to the minimally allowed value: 
by increasing the MAGN contribution from its minimal 
to 70% of its maximally expected value [27|, limits on a 
possible DM contribution could be improved by at least 
one order of magnitude compared to the most conser- 
vative assumptions concerning the astrophysical contri- 
butions (note that setting the 7-ray flux from MAGN 
to the maximally allowed value would actually overshoot 
the Fermi data by more than 2a). 

The right panel of Fig.[3l on the other hand, demon- 
strates that the uncertainty in the flux from other as- 
trophysical backgrounds generally has a less significant 
impact on the final limits. Here, we fix the MAGN con- 
tribution to its best-fit value [27[ and vary Ebms from 
its minimally to maximally predicted contribution as de- 
scribed in Section [TTJ We also show that the choice of the 
underlying SF galaxy model introduces a certain system- 



atic uncertainty: depending on the WIMP mass, the re- 
sulting limits arc up to a factor of ^3 stronger for the PL 
than for the MW model (unless Ebms only contributes 
at the minimal level, in which case there is no difference 
because the MAGN contribution always dominates). 

In Fig. m we show the effect of the MAGN contri- 
bution for various DM annihilation channels. In the 
top left panel, we present our most conservative limits. 
For bb, and t + t~ final states, and DM masses 

below ^100 GcV, they are comparable with the limits 
previously derived in Ref. (l7j . where a similar model- 
ing of the astrophysical background was considered; for 
m x > 100 GeV the inclusion of ICS becomes important 
and consequently significantly improves our new limits 
with respect to those for leptonic channels. The right 
panel shows a more realistic scenario, where we take 
all astrophysical backgrounds to contribute at their ex- 
pected 'best-fit' values. Compared to the conservative 
case, this shifts down the upper bounds by a factor of 
roughly 2. 

In the lower panels of Fig. @J we take a more opti- 
mistic approach in order to illustrate the full potential of 
using the IGRB to constrain galactic DM annihilation. 
Assuming MAGN to contribute at 60% of the maximal 
flux, for example, would result in limits that clearly start 
to approach the 'thermal' cross section of (av) ~ 3-10" 26 
cm 3 s _1 . Moving to the possibility of even higher levels 
for the MAGN contribution, we note that one quickly 
reaches a situation where already a small increase in the 
MAGN flux has an enormous impact on the allowed DM 
contribution - simply because the astrophysical back- 
grounds are already extremely close to the IGRB data 
points (with 2a errors). In this limit, the validity of 
our method to obtain upper bounds clearly breaks down 
and a more sophisticated statistical treatment would be 
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Conservative (min £ BS m and MAGN) Best fit (best fit I BSM and MAGN) 




M x [GeV] M z [GeV] 

FIG. 4: Upper limits on the velocity averaged DM annihilation cross section (av) into 7-rays from the Fermi-LAT IGRB data 
[TEj . for various possible annihilation channels. The top left (right) panel shows our most conservative (more realistic) limits, 
obtained by setting all astrophysical contributions to the minimally allowed ('best-fit') level. The bottom panels, where we 
also fix Ebms at the 'best-fit' value, represent a more optimistic choice for the MAGN contribution concerning potential DM 
limits: in the left (right) panel MAGN are assumed to actually contribute 60% (85%) of the maximally expected flux. Note 
that for the limits in the bottom right panel we require the total 7-ray flux not to overshoot any data point by 3a (rather than 
2a like in all other cases presented in this work). The grey area indicates the typical benchmark value for the cross section of 
thermally produced WIMPs, (av) ~ 3 • 1CT 26 cm 3 s -1 . 



needed. In order to illustrate that one may still poten- 
tially be able to constrain the 'thermal' annihilation rate 
up to DM masses of several TeV, we show in the bottom 
right panel of Fig. [3] the limits that result for an assumed 
MAGN contribution at 85% of the maximally predicted 
flux. Note that those limits were derived by requiring 
that no data point is exceeded by more than 3<r, rather 
than 2cr, in order to (partially) meet the above concerns. 

Let us stress in passing that the inclusion of ICS 7- 
rays is particularly relevant for large DM masses (m x > 
1 TeV) . This is the mass range which - for leptonic fi- 
nal states [72U75I - could in principle contribute to the 
observed rise in the cosmic ray positron fraction [76l - l78| . 
Our results strongly disfavor a DM candidate annihilat- 
ing mainly into e + e~ with (av) > 0.3 (1)-10 -23 cm 3 s _1 
for masses of about 1-3 (10) TeV. A bit looser constraints 
are obtained for annihilation into muons or tau leptons, 
the latter roughly at the level which is required to fit 



the AMS-02 data [zl [75|. Assuming a relatively large 
MAGN contribution, like e.g. in the scenario displayed 
in the bottom-left panel of Fig. [H the IGRB would thus 
allow to firmly rule out a DM explanation of the rising 
positron fraction. 

As already mentioned, our results improve the limits 
derived in Ref. [TtJ at DM masses above ^100 GeV due 
to the inclusion of ICS; the 'best-fit' scenario displayed 
in Fig. @] strengthens those limits further by a factor of 
^2 due to the contribution from MAGN. Our constraints 
for high masses are significantly stronger than those ob- 
tained from the ICS contribution alone in an analysis of 
preliminary Fermi data [451 ]. and stronger by an overall 
factor of roughly 4 (for bb) compared to those obtained 
for the analysis of ly IGRB data in a smaller sky region 
than considered here [79j . Compared to the Fermi anal- 
ysis of possible galactic DM contributions to the IGRB 
|3| , with no background modeling, our conservative limits 
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arc comparable for bb and t + t~ final states (for m x > 20 
GeV). The adopted background modeling in that analy- 
sis, on the other hand, leads to bounds comparable with 
our MAGN best fit choice for m x ~ 100 GeV; for higher 
masses m x > 1 TeV, the bounds in the MAGN best fit 
case are a factor of roughly 3 stronger. 

Let us finally remark that the bounds on (av) we have 
derived here are consistent with constraints from differ- 
ent targets in the sky obtained with very different as- 
trophysical backgrounds and measurement systematics. 
In particular, our 'optimistic' limits on the DM annihi- 
lation cross section - assuming an MAGN flux which is 
60% of its maximal theoretical prediction - are generally 
competitive with the bounds from the joint Likelihood 
analysis of 10 dwarf Spheroidal galaxies Q. In fact, for 
TeV masses and leptonic final states, our constraints are 
at least comparable even under the conservative assump- 
tion of a minimal MAGN contribution. We also obtain 
comparable results at low masses with bounds on lep- 
tonic channels from radio emission j80l] , which in turn 
depend significantly on astrophysical and cosmological 
assumptions. At masses m x > 100 GeV our bounds 
from 7-rays are stronger, independently of the amount of 
MAGN included in the background. Finally, our conser- 
vative constraints are stronger than those obtained from 
galaxy cluster Q and CMB observations [8l| . 

V. CONCLUSIONS 

In this work, we have considered the contribution of 
annihilating galactic DM to the IGRB and presented an 
updated calculation of the resulting upper bounds on the 
DM annihilation cross section (av). In deriving these 
bounds, we took into account both prompt and ICS emis- 
sion for DM- induced 7-rays. For the astrophysical 7- 
ray background, we included all relevant contributions 
at their minimally expected level, in particular the con- 
tribution from unresolved MAGN that has been pointed 
out only very recently. Given that we only consider high 
galactic latitudes, our results are very robust and barely 
depend on the choice of DM density profile (apart from 
an overall normalization corresponding to the local DM 
density, which may either strengthen or weaken those 
bounds by a factor of at most 3). Let us stress that our 
bounds are very conservative in that they do not include 
the expected - and potentially very sizable - effect from 
DM substructures at galactic or cxtragalactic distances. 

The resulting IGRB bounds on (av) are still competi- 
tive with, or more stringent than, limits from radio data, 
CMB or galaxy cluster observations. In particular, our 
results - along with stringent limits from antiproton data 
[12-1111 - strongly disfavor a DM-based explanation of 
the unexpected rise in the cosmic ray positron fraction 



at high energies that was recently confirmed by AMS-02 
[7§ | , at least if MAGN do not contribute at the most min- 
imal level currently allowed. The main reason for this is 
the sizable ICS emission in 7-rays that is induced by the 
high annihilation rate into leptonic channels required to 
fit the positron data. 

Even more importantly, we have demonstrated that 
there is significant room for future improvement of limits 
on the DM annihilation rate derived from the IGRG. In 
fact, under the optimistic (though by no means unrealis- 
tic) assumption that the astrophysical processes consid- 
ered here contribute more or less at the maximal level 
consistent with our present understanding, a 'thermal' 
value of 3 • 10~ 26 cm 3 s _1 for the annihilation cross sec- 
tion could be excluded up to DM masses of several TeV. 
Notably, such a limit will be impossible to reach even 
for far-future ground-based observations of point- or ex- 
tended sources like dwarf galaxies or the galactic cen- 
ter. Even for space-based observations of dwarf galax- 
ies, which presently sets the most robust and stringent 
bounds of this kind, it will be extremely challenging for 
extended future missions to probe thermally produced 
DM for such large masses. 

The most important contribution to the IGRB origi- 
nates from unresolved MAGN, which we have included 
here for the first time in the context of putting limits 
on DM annihilation. In order to actually exploit the 
great potential of the IGRB discussed above, a reduction 
of the relatively large theoretical uncertainty on the dif- 
fuse 7-ray flux from unresolved MAGN is thus manda- 
tory. In particular, one would need several tens of de- 
tected MAGN in 7-rays to be able to directly determine 
their 7-ray luminosity function (which is presently de- 
rived from radio observations); this goal is well in reach 
for Fermi-hAT and next generation 7-ray telescopes. At 
some point, also an improved understanding of the other 
astrophysical contributions will help to push limits fur- 
ther down - especially if one can develop the underlying 
modeling to a level that allows to predict those contri- 
butions not only as an all-sky average but as a function 
of the viewing angle. Given the expected significantly 
increased statistics for future IGRB observations, finally, 
chances are good that the IGRB will turn out to be one 
of the most efficient future means of probing thermally 
produced DM. 

Acknowledgments 

We warmly thank A. Bottino, N.Fornengo and L. La- 
tronico for useful discussions and a careful reading of 
the manuscript. T.B. and F.C. gratefully acknowledge 
support from the German Research Foundation (DFG) 
through Emmy Noether grant BR 3954/1-1. 



[1] T. Bringmann and C. Weniger, Phys.Dark Univ. 1, 194 [2] W. Atwood et al. (LAT Collaboration), Astrophys.J. 
(2012), 1208.5481. 



9 



697, 1071 (2009), 0902.1089. [31 

M. Ackermann, M. Ajello, W. B. Atwood, L. Baldini, 

G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, [32 

R. D. Blandford, E. D. Bloom, et al., Astrophys. J. 761, 

91 (2012), 1205.6474. 

A. A. Abdo, M. Ackermann, M. Ajello, W. B. Atwood, 
L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, K. Bech- [33 
tol, R. Bellazzini, et al., Astrophys. J. 712, 147 (2010), 
1001.4531. 

M. Ackermann, M. Ajello, A. Allafort, L. Baldini, J. Bal- 
let, G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, [34 
R. D. Blandford, et al., JCAP 5, 25 (2010), 1002.2239. 

A. A. Abdo, M. Ackermann, M. Ajello, L. Baldini, J. Bal- 
let, G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, 

B. Berenji, et al., JCAP 4, 14 (2010), 1002.4415. [35 
M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, 

L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, K. Bech- 
tol, R. Bellazzini, et al., Phys. Rev. Lett. 107, 241302 

(2011) , 1108.3546. [36 
T. Bringmann, X. Huang, A. Ibarra, S. Vogl, and 

C. Weniger, JCAP 7, 54 (2012), 1203.1312. [37 

C. Weniger, JCAP 8, 7 (2012), 1204.2797. 

E. Tempel, A. Hektor, and M. Raidal, JCAP 1209, 032 [38 

(2012) , 1205.1045. 

D. P. Finkbeiner, M. Su, and C. Weniger, JCAP 1, 29 [39 

(2013) , 1209.4562. 

E. Bloom, E. Charles, E. Izaguirre, A. Snyder, A. Albert, 
et al. (2013), 1303.2733. 

Fermi-LAT Collaboration, ArXiv e-prints (2013), [40 
1305.5597. 

C. Weniger, M. Su, D. P. Finkbeiner, T. Bringmann, and [41 
N. Mirabal (2013), 1305.4710. 

A. A. Abdo and et al. et al., Phys. Rev. Lett. 104, 101101 [42 
(2010), 1002.3603. 

M. Ackermann et al. (Fermi LAT Collaboration), [43 
Phys.Rev. D85, 083007 (2012), 1202.2856. 

F. Calore, V. de Romeri, and F. Donato, Phys. Rev. D [44 
85, 023004 (2012), 1105.4230. 
L. Bergstrom, J. Edsjo, and P. Ullio, Phys. Rev. Lett. 87, [45 
251301 (2001), astro-ph/0105048. 

P. Ullio, L. Bergstrom, J. Edsjo, and C. G. Lacey, [46 
Phys.Rev. D66, 123502 (2002), astro-ph/0207125. 
J. E. Taylor and J. Silk, Mon.Not.Roy.Astron.Soc. 339, [47 
505 (2003), astro-ph/0207299. 

M. Papucci and A. Strumia, JCAP 1003, 014 (2010), [ 
0912.0742. [49 
M. Cirelli, P. Panci, and P. D. Serpico, Nucl.Phys. B840, 
284 (2010), 0912.0663. [50 

E. J. Baxter, S. Dodelson, S. M. Koushiappas, and L. E. 
Strigari, Phys.Rev. D82, 123511 (2010), 1006.2399. 

K. N. Abazajian, S. Blanchet, and J. P. Harding, [51 
Phys.Rev. D85, 043509 (2012), 1011.5090. 

S. Blanchet and J. Lavalle, JCAP 1211, 021 (2012), [52 
1207.2476. 

J. Zavala, V. Springel, and M. Boylan-Kolchin, 
Mon.Not.Roy.Astron.Soc. 405, 593 (2010), 0908.2428. 

F. Calore, M. Di Mauro, F. Donato, M. Ajello, and L. La- [53 
tronico, Submitted to Astrophys. J. (2013). 
M. Di Mauro, Talk at Moriond 2013, Very High Energy [54 
Phenomena in the Universe. La Thuile, March 10 (2013). 

A. Abramowski et al. (H.E.S.S. Collaboration), [55 
Phys.Rev.Lett. 106, 161301 (2011), 1103.3266. 

B. Acharya, M. Actis, T. Aghajani, G. Agnetta, [56 
J. Aguilar, et al., Astropart.Phys. 43, 3 (2013). 



C. M. Urry and P. Padovani, PASP 107, 803 (1995), 
arXiv:astro-ph/9506063. 

A. A. Abdo, M. Ackermann, M. Ajello, A. Allafort, 
E. Antolini, W. B. Atwood, M. Axelsson, L. Baldini, 
J. Ballet, G. Barbiellini, et al., Astrophys. J. S.S. 188, 
405 (2010), 1002.2280. 

P. L. Nolan, A. A. Abdo, M. Ackermann, M. Ajello, 
A. Allafort, E. Antolini, W. B. Atwood, M. Axelsson, 
L. Baldini, J. Ballet, et al., Astrophys. J. S.S. 199, 31 
(2012), 1108.1435. 

A. A. Abdo, M. Ackermann, M. Ajello, E. Antolini, 
L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, B. M. 
Baughman, K. Bechtol, et al., Astrophys. J. 720, 435 

(2010) , 1003.0895. 

M. Ajello, M. S. Shaw, R. W. Romani, C. D. Dermer, 
L. Costamante, O. G. King, W. Max-Moerbeck, A. Read- 
head, A. Reimer, J. L. Richards, et al., Astrophys. J. 
751, 108 (2012), 1110.3787. 

A. K. Harding, V. V. Usov, and A. G. Muslimov, Astro- 
phys. J. 622, 531 (2005), arXiv:astro-ph/0411805. 
C.-A. Faucher-Giguere and A. Loeb, JCAP 1, 5 (2010), 
0904.3102. 

The Fermi-LAT collaboration, ArXiv e-prints (2013), 
1305.4385. 

M. Ackermann, M. Ajello, A. Allafort, L. Baldini, J. Bal- 
let, D. Bastieri, K. Bechtol, R. Bellazzini, B. Berenji, 

E. D. Bloom, et al., Astrophys. J. 755, 164 (2012), 
1206.1346. 

F. W. Stecker and T. M. Venters, Astrophys. J. 736, 40 

(2011) , 1012.3678. 

H. Bengtsson, P. Salati, and J. Silk, Nuclear Phys. B 
346, 129 (1990). 

L. Bergstrom, P. Ullio, and J. H. Buckley, Astropart. 
Phys. 9, 137 (1998), 9712318. 

A. Bottino, F. Donato, N. Fornengo, and S. Scopel, Phys. 

Rev. D 70, 015005 (2004), 0401186. 

T. Sjostrand, S. Mrenna, and P. Z. Skands, Com- 

put.Phys.Commun. 178, 852 (2008), 0710.3820. 

M. Cirelli and P. Panci, Nuclear Physics B 821, 399 

(2009), 0904.3830. 

M. Regis and P. Ullio, Phys.Rev. D78, 043505 (2008), 
0802.0234. 

X. H. Sun, W. Reich, A. Waelkens, and T. A. Enfilin, 

A&A 477, 573 (2008), 0711.1572. 

J. L. Han and G. J. Qiao, A&A 288, 759 (1994). 

M. Prouza and R. Smi'da, A&A 410, 1 (2003), 

arXiv:astro-ph/0307165. 

M. S. Pshirkov, P. G. Tinyakov, P. P. Kronberg, and 
K. J. Newton-McGee, Astrophys. J. 738, 192 (2011), 
1103.0814. 

A. Burkert, Astrophys. J. Lett. 447, L25 (1995), 
arXiv:astro-ph/9504041. 

A. J. Deason, V. Belokurov, N. W. Evans, S. E. Ko- 
posov, R. J. Cooke, J. Peharrubia, C. F. P. Laporte, 
M. Fellhauer, M. G. Walker, and E. W. Olszewski, 
Mon.Not.R.Astron.Soc. 425, 2840 (2012), 1205.6203. 
O. Y. Gnedin, W. R. Brown, M. J. Geller, and S. J. 
Kenyon, Astrophys. J. 720, L108 (2010), 1005.2619. 
A. Irrgang, B. Wilcox, E. Tucker, and L. Schiefelbein, 
A&A 549, A137 (2013), 1211.4353. 

P. R. Kane, S. Sharma, G. F. Lewis, and J. Bland- 
Hawthorn, Astrophys. J. 761, 98 (2012), 1210.7527. 
C. S. Kochanek, Astrophys. J. 457, 228 (1996), 
arXiv:astro-ph/9505068. 



10 



[57] T. Sakamoto, M. Chiba, and T. C. Beers, A&A 397, 899 
(2003), arXiv:astro-ph/0210508. 

[58] L. L. Watkins, N. W. Evans, V. Belokurov, M. C. 
Smith, P. C. Hewett, D. M. Bramich, G. F. Gilmore, 
M. J. Irwin, S. Vidrih, L. Wyrzykowski, et al., 
Mon.Not.R.Astron.Soc. 398, 1757 (2009), 0906.0498. 

[59] M. I. Wilkinson and N. W. Evans, 
Mon.Not.R.Astron.Soc. 310, 645 (1999), arXiviastro- 
ph/9906197. 

[60] X. X. Xue, H. W. Rix, G. Zhao, P. Re Fiorentin, T. Naab, 
M. Steinmetz, F. C. van den Bosch, T. C. Beers, Y. S. 
Lee, E. F. Bell, et al., Astrophys. J. 684, 1143 (2008), 
0801.1232. 

[61] J. F. Navarro, C. S. Frenk, and S. D. White, Astrophys. J. 

462, 563 (1996), astro-ph/9508025. 
[62] J. Einasto, Trudy Astrofizicheskogo Instituta Alma-Ata 

5, 87 (1965). 

[63] J. I. Read, G. Lake, O. Agertz, and V. P. Debattista, 

Mon. Not. R. Astron. Soc. 389, 1041 (2008), 0803.2714. 
[64] J. Diemand and B. Moore, Advanced Science Letters 4, 

297 (2011), 0906.4340. 
[65] V. Springel, J. Wang, M. Vogelsberger, A. Ludlow, 

A. Jenkins, A. Helmi, J. F. Navarro, C. S. Frenk, 

and S. D. M. White, Mon.Not.R.Astron.Soc. 391, 1685 

(2008), 0809.0898. 
[66] L. Bergstrom, J. Edsjo, and P. Ullio, Phys. Rev. Lett. 

87, 251301 (2001), arXiv:astro-ph/0105048. 
[67] P. Ullio, L. Bergstrom, J. Edsjo, and C. Lacey, Phys. 

Rev. D 66, 123502 (2002), arXiv:astro-ph/0207125. 
[68] J. Diemand, M. Kuhlen, and P. Madau, Astrophys. J. 

657, 262 (2007), astro-ph/0611370. 



[69] M. Kuhlen, J. Diemand, and P. Madau, Astrophys. J. 

686, 262 (2008), 0805.4416. 
[70] L. Pieri, J. Lavalle, G. Bertone, and E. Branchini, 

Phys.Rev. D83, 023518 (2011), 0908.0195. 
[71] P. Salucci, F. Nesti, G. Gentile, and C. Frigerio Martins, 

A&A 523, A83 (2010), 1003.3101. 
[72] L. Bergstrom, J. Edsjo, and G. Zaharijas, Phys. Rev. Lett. 

103, 031103 (2009), 0905.0333. 
[73] Q. Yuan, X.-J. Bi, G.-M. Chen, Y.-Q. Guo, S.-J. Lin, 

et al. (2013), 1304.1482. 
[74] I. Cholis and D. Hooper (2013), 1304.1840. 
[75] H.-B. Jin, Y.-L. Wu, and Y.-F. Zhou (2013), 1304.1997. 
[76] O. Adriani et al. (PAMELA Collaboration), Nature 458, 

607 (2009), 0810.4995. 
[77] M. Ackermann et al. (Fermi LAT Collaboration), 

Phys.Rev.Lett. 108, 011103 (2012), 1109.0521. 
[78] M. Aguilar et al. (AMS Collaboration), Phys.Rev.Lett. 

110, 141102 (2013). 
[79] M. Cirelli, P. Panci, and P. D. Serpico, Nuclear Physics 

B 840, 284 (2010), 0912.0663. 
[80] N. Fornengo, R. A. Lineros, M. Regis, and M. Taoso, 

JCAP 1, 5 (2012), 1110.4337. 
[81] S. Galli, F. Iocco, G. Bertone, and A. Melchiorri, Phys. 

Rev. D 84, 027302 (2011), 1106.1528. 
[82] F. Donato, D. Maurin, P. Brun, T. Delahaye, and 

P. Salati, Phys. Rev. Lett. 102, 071301 (2009), 

0810.5292. 

[83] M. Cirelli and G. Giesen, JCAP 4, 015 (2013), 1301.7079. 
[84] C. Evoli, I. Cholis, D. Grasso, L. Maccione, and P. Ullio, 
Phys. Rev. D 85, 123511 (2012). 



